Spin nematics in the bilinear-biquadratic S=l spin chain 
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We report the existence of an extended critical, nondimerized region in the phase diagram of 
the bilinear-biquadratic spin-one chain. The dominant power law correlations are ferroquadrupolar 
i.e. spin nematic in character. Another known critical region is also characterized by dominant 
quadrupolar correlations, although with a different wave vector. Our results show that spin ne- 
matic correlations play an important role in quantum magnets with spin S > 1 in regions between 
antiferromagnetic and ferromagnetic phases. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Mg, 75.40.Cx 



The recent experimental demonstration of the tran- 
sition to a Mott insulating state of atoms in an opti- 
cal lattice from a superfluid state has opened the way 
to novel realizations of effective quantum lattice mod- 
els with widely tunable control parameters. Quantum 
magnetic systems can be realized by spinor atoms in an 
optical lattice, e.g. ^'^Na with a total S = 1 moment. 
Confining S = \ atoms to an optical lattice there are two 
scattering channels for identical atoms with total spin 
S = Q,2 which can be mapped to an effective bilinear 
and biquadratic spin interaction [^|^: 



H 



Jq (Si • Sj) 



(1) 



where we adopt the standard parametrization J = 
cos 9 and Jq — svnO. In one dimension, the bilinear- 
biquadratic spin-one model has a rich phase diagram (see 
Fig.GJ, which includes the well-known Haldane gap phase 
^] , a dimerized phase , a ferromagnetic phase and - as 
we will show in this Letter - two critical spin nematic 
phases. In a pioneering paper |Q Chubukov suggested 
the existence of a gapped, nondimerized phase with dom- 
inant spin nematic correlations close to the ferromagnetic 
region of the phase diagram. Subsequent numerical work 
13 could however not substantiate this claim and it was 
therefore believed that the dimerized phase would ex- 
tend up to the ferromagnetic phase boundary. Recent 
quantum Monte Carlo calculations and field theoret- 
ical work ^] suggest that this picture might need to be 
reconsidered, especially in the light of possible experi- 
mental verifications in Bose-Einstein condensate systems 

. . . 

In dimensions higher than one, the bilinear-biquadratic 

spin-one model is well understood in the relevant param- 
eter regime. It exhibits a spin nematic phase with a finite 
spin quadrupole moment in the groundstate 0, 0| • 

In this Letter we propose a novel scenario for the phase 
diagram of the one-dimensional model in between the 
dimerized and the ferromagnetic phases. In particular, 
we substantiate the existence of an extended critical, 



nondimerized region with dominant spin nematic corre- 
lations using extensive numerical simulations. We fur- 
ther characterize a second extended critical phase be- 
tween the Haldane and the ferromagnetic phase, as also 
having prevalent quadrupolar correlations, albeit with a 
different wavevector. 

The outline of this paper is as follows. In the first 
part we present various numerical results providing ev- 
idence for an extended critical phase without sponta- 
neous dimerization for 6 G (— 37r/4, — O.BTtt]. We show 
that the quintuplet gap vanishes over this whole interval 
and that the transition between the dimerized and the 
critical phase is of a generalized Beresinski-Kosterlitz- 
Thouless (BKT) type. The absence of finite dimerization 
is discussed. We further determine conformal field the- 
ory parameters, such as scaling dimensions and the cen- 
tral charge for the critical region. Finally we address the 
second extended critical phase 9 G [+-K /A,+-k /2), which 
we exhibit to have prevailing quadrupolar correlations at 
wavevectors ±27r/3. 

Numerical techniques — Our results constitute the 



critical, 

/c=±27t/3 quadrupolar 




critical, 
k=Q quadrupolar 



FIG. 1: Phase diagram of the bilinear-biquadratic spin-one 
chain. The new elements are (a) an extended, critical phase 
with dominant ferroquadrupolar correlations for — 37r/4 < 
9 < — O.GTtt and (b) the characterization of the other criti- 
cal region 7r/4 < 6^ < 7r/2 as having dominant k — ±27r/3 
quadrupolar correlations. 
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FIG. 2; Extrapolated DMRG S = 2 gaps as a function of 
9. The gray region denotes the estimated onset of the critical 
behavior. The filled diamond is the exactly known gap result 
a.t 6 — — 7r/2. Inset; finite size extrapolation of the S — 2 
gaps (32 to 512 sites). 9/tt varies from —0.5 to —0.7 with 
decrements of 0.05 from top to bottom. 

synthesis of different computational methods: i) Exact 
diagonalization (ED) of periodic chains up to 18 sites, 
using both static and dynamic observables, ii) Large 
scale Density Matrix Renormalization Group (DMRG) 
[l^ calculations using the finite size algorithm on systems 
up to 512 sites while keeping 1000 states, iii) Strong cou- 
pling series expansions 13] around the dimerized limit 
up to 10th order. 

Closing of the Gap — We first track the evolution of 
the energy gap to the first excitation using DMRG cal- 
culations on long open chains. A phase transition is sig- 
naled by the closing of the gap. The lowest excited state 
on open chains carries spin two for 9 G (— 37r/4, —tt/2] 
0. Our results shown in Fig. [21 are clearly consistent 
with a finite gap of the S = 2 excitation for 9 > — 0.657r. 
However, for 9 < — O.GStt the extrapolated gap is very 
small (of the order of 10""^ for 9 — — O.Ttt) which sug- 
gests the existence of a phase transition around — O.BTtt 
below which the gap is zero. The data shows no evi- 
dence for a reopening of the S = 2 gap in the interval 
9 G (— 37r/4, — 0.677r], at variance with the initial pro- 
posal . In order to investigate the closing of the gap by 
different means we have calculated strong coupling series 
expansions of the S" = 2 single particle gap starting from 
the dimerized limit up to 10th order in the interdimcr 
coupling A for fixed 9. A direct evaluation of the quintu- 
plet gap shows a closing of the gap around 9 = — O.GTtt 
which is illustrated in the left panel of Fig. |3| Using 
DLog Fade approximants we have determined the crit- 
ical interdimer coupling Xc{9) where the gap closes, see 
upper right panel of Fig. |3| Consistent with previous es- 
timates by Chubukov [6j we find that the gap closes for 
the uniformly coupled (A = 1) chain around 9 « — O.GTtt. 

Having substantiated the closure of the gap in both 
DMRG and series expansions with indications of an ex- 
tended gapless region — 37r/4 < 9 < — O.GTtt, we now con- 
sider the existence of a generalized BKT phase transition. 
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FIG. 3: Strong coupling dimer expansion of the S = 2 gap. 
Left panel: 5 = 2 gap as a function of 9. The gap closes 
at 9c ~ — 0.677r. Upper right panel: Critical value Ac where 
the gap vanishes. The uniformly coupled chain of interest 
here corresponds to A = 1 (dashed line). Lower right panel: 
Critical exponent of the S — 2 gap opening as a function of 
the dimerization (A ~ |A — Xcp) calculated by various DLog 
Fade approximants. 
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FIG. 4: Finite size scaling of the crossing points between 
the singlet gap at momentum n and the quintuplet {S = 2) 
gap at momentum 0. The extrapolation yields a critical value 
9c ~ — 0.677r in the thermodynamic limit. 

To determine the critical point we use phenomenological 
level spectroscopy. We calculate the level crossing 9^ of 
the lowest singlet excitation at fc = tt with the lowest 
spin-2 level at A: = for different system sizes L with ED 
and extrapolate L — > oo. The results are shown in Fig.^ 
The extrapolation is performed with 1/L and cor- 
rections, fitting only the largest four system sizes. The 
extrapolated critical point is 9c = (—0.67 ± 0.02)7r, in 
good agreement with estimates obtained by other tech- 
niques 0, 01 . The coefficient of the 1/L correction is 
small, and fitting without this term yields slightly larger 
values 9c > — O.GTtt. 

Dimerization and quadrupolar correlations — Dimer- 
ization in the phase diagram of the spin-one chain has 
been firmly established by exact results obtained for 
the 9 = — 7r/2 point j^. It should be noted that al- 
though the system is dimerized, it is rather poorly de- 
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FIG. 5: Static and generalized susceptibilities for the bilin- 
ear dimerization (left panel) and the ferroquadrupolar corre- 
lations (right panel), obtained by ED on systems of 8 to 16 
sites. Solid symbols: Static structure factors. Open symbols: 
Corresponding generalized susceptibilities. 

scribed by a simple product wavefunction of alternat- 
ing singlets, as can be seen from the small gap and the 
large correlation length over the whole dimerized phase. 
The dimerization order parameter |l4l | considered here 
is: {{Si-i ■ Si) — (Si ■ Si+i)). The second kind of cor- 
relations expected to be important in this region of the 
phase diagram [lol are the quadrupolar spin fluctua- 
tions C^Ur) ^ {[{S^r " 2/3] [(5,^)^ - 2/3]). 

We have calculated the standard static structure factor 
(SF) and the generahzed nonhnear susceptibility (GNS) 
|l5| for both kinds of correlations. We present ED 
results regarding the dimerization in the left panel of 
Fig. □ The SF and the GNS both diverge as L ^ oo 
for 9 = — 7r/2, where long range dimer order is estab- 
lished Similar behavior is found for > — 0.677r. For 
9 e (-37r/4, -0.677r) however, we find that the SF and 
the GNS are both decreasing with system sizes, indicating 
the absence of static dimerization in this gapless region. 
In the right panel we display the corresponding quantities 
for the k — (ferroquadrupolar) mode of the quadrupolar 
correlations. Here the behavior is different: while there 
are only short range correlations deep in the dimerized 
phase, the ferroquadrupolar correlations increase drasti- 
cally with system sizes for 6 £ (— 37r/4, — O.GTtt). Since 
quadrupolar long ranged order is forbidden at T = in 
this ID model - apart from 9 — — 37r/4, -|-7r/4, where 
the order parameter commutes with the Hamiltonian - 
we therefore expect a power-law behavior. Note that 
exactly at 9 = — 37r/4 the system is a SU{3) fcrromag- 
net and shows ferromagnetism and ferroquadrupolar long 
range order at the same time [l^ . 

Finally we want to discuss the relation of the proposed 
scenario to the numerical results presented in Ref. 
where an intermediate phase as proposed by Chubukov 
could not be found. We agree with the authors on the 
absence of a gapped, nematic-like phase. However, we 
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FIG. 6: Conformal field theory parameters calculated in ED 
on chains of 8 up to 18 sites. The scaling dimensions of the 
singlet (fc = tt), and the quintuplet [k — 0) field are shown. 
Inset: effective central charge. The value at the boundary of 
the critical region (6 ~ — 0.677r) is c « 3/2. 

think that the numerical results in 7] are compatible with 
the existence of an intermediate critical, nondimerized 
phase, though the authors have not explicitly considered 
such a phase. 

Field theoretical description — In order to discrimi- 
nate between potential conformal theories describing the 
phase transition and the extended gapless region, we 
calculate relevant field theory parameters, i.e. the scal- 
ing dimensions (SD) x of several fields and the central 
charge c (Fig. The calculation of these parameters 
has been performed with ED and relies on finite size scal- 
^ properties of groundstate and excited states energies 
Leading logarithmic corrections have been taken 
into account. The SD of the k = tt singlet field {xg) and 
the fc = quintuplet field (xq) are both close to 3/8 at 
the onset of criticality near 9 « — 0.677r. Xg significantly 
increases, while Xq decreases and seems to approach 
a.s 9 ^ —371-/4, compatible with quadrupolar long range 
order precisely at 9 = — 37r/4 Our results for the ef- 
fective central charge are unexpected. The central charge 
c does not seem to be constant throughout the critical re- 
gion. The smallest value is found at the onset (c « 3/2), 
and then seems to increase monotonously. While this 
behavior can not directly be ruled out on field theoret- 
ical grounds, it is rather uncommon. It remains to be 
understood whether c is really continuously increasing 
or whether we are facing a crossover phenomenon. The 
critical theory at 9 fa — 0.677r we however believe to be 
well characterized by a level two SU{2) Wess-Zumino- 
Witten model given both the scahng dimensions (3/8) 
and the central charge (3/2). Further support for this 
claim comes from the indirect calculation of Xg within 
our series expansions. There the critical exponent 7 of 
the 5 = 2 gap is related to the scaling dimension Xs by 
7 = 1/(2 — Xs)- In the lower right panel of Fig. 13 we 
show the critical exponent calculated by DLog Pade ap- 
proximants to the 10th order series. In the vicinity of 
9 w — 0.677r the various approximants reveal only small 
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FIG. 7: Static structure factors for spin and quadrupolar 
correlations at momentum fc = 27r/3 as a function of 6 for 
different system sizes. At the SU(3) point 9 = +7r/4 tlie two 
correlation functions are symmetry related. By going deeper 
into the critical region the spin correlations get weaker, while 
the quadrupolar correlations are enhanced. Inset: the central 
charge, the excitation velocity and the scaling dimension of 
the quintuplet field in the critical region. 

spreading. Within the precision of a 10th order calcula- 
tion the exponents comply with a scaling dimension of 
Xs — 3/8 (dashed line). 

The period 3 phase — The existence of an extended 
critical phase with soft modes at A: = ±27r/3 has first 
been reported in llSl and further analyzed and inter- 
preted in 0, 123, 12 ij . Our results presented in the first 
part raise the question whether this critical region is of 
spin nematic character as well. Further hints come from 
the field theoretical work of Itoi and Kato j^^, which 
shows that throughout the critical region the spin 2 mode 
at A: = ±27r/3 has correlation-enhancing logarithmic cor- 
rections. This is consistent with the earlier observation, 
that the lowest level at k = ±27r/3 carries spin 2 We 
have calculated the static spin structure factor [S = 1) 
and the spin quadrupolar structure factor {S = 2) both 
at momentum 27r/3. The results shown in Fig. |7| display 
nicely that the quadrupolar correlations become more im- 
portant than the corresponding spin-spin correlations as 
one leaves the SU{i) symmetric point. Although all cor- 
relation exponents (dimer, spin, quadrupolar) are equal 
to 77 = 4/3, logarithmic corrections render the quadrupo- 
lar correlations dominant in this phase. Note that this 
phenomenon is similar to the 5=1/2 antiferromagnetic 
Heisenberg chain, where both spin-spin and dimer-dimer 
correlations have equal correlation exponents rj — 1^ but 
the logarithmic corrections choose the spin-spin correla- 
tions to decay more slowly. For completeness we show 
selected field theory parameters for this critical region 
in the inset of Fig. The results are in good agreement 
with field theoretical predictions [2l'| and support the ap- 
plicability of the used ED techniques to determine field 
theory parameters. 

Conclusions and Outlook — To summarize, we have 
shown that spin nematic correlations play an impor- 



tant role in the phase diagram of the S = 1 bilinear- 
biquadratic chain. Applying a combination of numerical 
techniques we reconsidered Chubukov's proposal of 
a gapped, nondimerized and nematic-like phase between 
the dimerized and the ferromagnetic phases. Our nu- 
merical results provide evidence for a modified picture 
with a phase transition at « — O.GTtt where the dimer- 
ized phase terminates. However, the phase in the win- 
dow E (— 37r/4, — 0.677r] remains critical, i.e. gapless 
and nondimerized. The dominant correlations are spin- 
quadrupolar correlations at A: = 0. In the second crit- 
ical region for 9 E [+7r/4, -f 7r/2) we also find dominant 
quadrupolar correlations at wavevectors k — ±27r/3. It 
will be interesting to compare our field theoretical pa- 
rameters to future analytical calculations. 
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